*******************************************************************************
* appendix figure 13: near-far all incomes rent event studies

********************************************************************************

cd "$main_dir/data"

* 250m event studies
local max_dist 600
local iso_dist 250
local g ng

use near_far_m`max_dist'_i`iso_dist'_`g',clear

egen ring_fe=group(building_id l250)
gen zero=0

reghdfe lrent tx_n3 tx_n2 zero tx_0 tx_1 tx_2 tx_3 /// 
	i.listing_bedroom i.listing_bathroom ///
	if building_sample==1 &  tx250_min_number>0 & tx256_min_number>0 ///
	[aweight=building_sample_weight] ,absorb(ring_fe building_id_year) vce(cl building_id)

coefplot, ///
	keep(tx* zero) ///			
	vert ///
	yline(0) xtitle(Years to treatment) ///
	rename(tx_n3 = "-3" tx_n2 = "-2" zero = "-1" tx_0="0" tx_1 = "1" tx_2 = "2" tx_3 = "3" tx_4="4")  ///
	ytitle("Event Study Coefficient") graphregion(color(white)) bgcolor(white) ///
	yline(0) omitted

* 400m event studies
local max_dist 600
local iso_dist 400
local g ng

use near_far_m`max_dist'_i`iso_dist'_`g',clear

egen ring_fe=group(building_id l400)
gen zero=0

reghdfe lrent tx_n3 tx_n2 zero tx_0 tx_1 tx_2 tx_3 /// 
	i.listing_bedroom i.listing_bathroom ///
	if building_sample==1 &  tx400_min_number>0 & tx46_min_number>0 ///
	[aweight=building_sample_weight] ,absorb(ring_fe building_id_year) vce(cl building_id)

coefplot, ///
	keep(tx* zero) ///			
	vert ///
	yline(0) xtitle(Years to treatment) ///
	rename(tx_n3 = "-3" tx_n2 = "-2" zero = "-1" tx_0="0" tx_1 = "1" tx_2 = "2" tx_3 = "3" tx_4="4")  ///
	ytitle("Event Study Coefficient") graphregion(color(white)) bgcolor(white) ///
	yline(0) omitted
